Pygmy dipole mode in deformed neutron-rich Mg isotopes close to the drip line 
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We investigate the microscopic structure of the low-lying isovector-dipole excitation mode in 
neutron-rich 36 ' 38,40 Mg close to the drip line by means of the deformed quasiparticle-random-phase- 
approximation employing the Skyrme and the local pairing energy-density functionals. It is found 
that the low-lying bump structure above the neutron emission-threshold energy develops when the 
drip line is approached, and that the isovector dipole strength at E x < 10 MeV exhausts about 6.0% 
of the classical Thomas-Reiche-Kuhn dipole sum rule in 40 Mg. We obtained the collective dipole 
modes at around 8 — 10 MeV in Mg isotopes, which consist of many two-quasiparticle excitations of 
neutron. The transition density clearly shows an oscillation of the neutron skin against the isoscalar 
core. We found significant coupling effects between the dipole and octupole excitation modes due 
to the nuclear deformation. It is also found that the responses for the compressional dipole and 
isoscalar octupole excitations are much enhanced in the lower energy region. 

PACS numbers: 21.10.Re; 21.60.Ev; 21.60.Jz 



I. INTRODUCTION 

Nuclei far from the stability has attracted a consider- 
able interest both experimentally and theoretically. Ex- 
ploring the multipole responses in unstable nuclei is in 
particular of great interest because they provide infor- 
mation on collective modes of excitation. In neutron- 
rich nuclei, the surface structure is quite different to the 
stable ones due to the presence of the loosely bound 
neutrons. One of the unique structures is the neutron 
skin [l], 0]. Since the collective excitations are sensi- 
tive to the surface structure, we can expect new kinds 
of exotic excitation mode associated with the neutron 
skin to appear in neutron-rich nuclei. An example is 
the soft dipole excitation or the pygmy dipole reso- 
nance (PDR), which is observed not only in light halo nu- 
clei 3, H, EH, H , H, [9l [Tol Ull . H2I Il3j . but also in heavier sys- 
tems [1J, [La [lg, LL2|> where an appreciable El strength 
is observed above the neutron emission threshold, and 
it exhausts several percents of the energy-weighted sum 
rule (EWSR) value. 

The structure of the PDR and its collectivity has 
been studied based on the mean-field calculations by 
many groups. (See Ref. [H| for extensive lists of ref- 
erences concerning the theoretical investigation of the 
PDR.) These studies, however, are largely restricted 
to spherical systems. Quite recently, by means of 
the deformed quasiparticle-random-phase approximation 
(QRPA), the low- lying dipole excitation in neutron-rich 
Ne isotopes [l9| and in neutron-rich Sn isotopes [2(| have 
been investigated. 

Presently, small excitation energies of the first 2+ 
state and striking enhancements of B(E2;0f — > 2^) in 
32 Mg [H HI and 34 Mg [H [H, [U are under lively 
discussions in connection with onset of the quadrupole 
deformation, breaking of the = 20 spherical magic 
number, pai ring co rrelation and continuum coupling ef- 
fects [H, [23, 128l.l29l. [30 | . In order to get clear understand- 
ing of the nature of quadrupole deformation and pairing 



correlations, it is strongly desirable to explore, both ex- 
perimentally and theoretically, excitation modes in Mg 
isotopes toward a drip line [111, [H, [H, [H, HH, [H, H3] • 

In the present article, we investigate the microscopic 
structure of the low-lying dipole excitation in neutron- 
rich Mg isotopes close to the drip line, paying special 
attention to the deformation effects on them. In a de- 
formed system, the soft negative-parity excitation modes 
could emerge associated with coupling between the dipole 
and octupole modes of excitation [381 ]. To this end, 
we perform a deformed QRPA in the matrix formal- 
ism on top of the coordinate-space Skyrme-Hartree-Fock- 
Bogoliubov (HFB) theory developed in Ref. 39]. The 
matrix formalism of the QRPA is suitable to investigate 
the microscopic structure of the excitation modes. 

This article is organized as follows: In the next section, 
the deformed Skyrme-HFB + QRPA method is recapit- 
ulated. In Sec. IIIIl we show the results of the deformed 
QRPA and investigate microscopic structures of the low- 
lying dipole states in 36 > 38 > 40 Mg. We discuss properties of 
the coupling among different modes of excitation unique 
in deformed neutron-rich nuclei in Sec. IIVI The last sec- 
tion is devoted to a summary. 



II. MODEL 

A detailed discussion of the deformed Skyrme- HFB + 
QRPA can be found in Ref. [3!|. Therefore, we just 
briefly recall the outline of the calculation scheme. 

In order to describe simultaneously the nuclear defor- 
mation and the pairing correlations including the un- 
boun d q uasiparticle states, we solve the HFB equa- 
tions do, S3 



h q - X q 
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in coordinate space using cylindrical coordinates r 
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(p,z,4>). We assume axial and reflection symmetries. 
Here, q = v (neutron) or ir (proton). For the mean- 
field Hamiltonian h, we employ the SkM* interaction [42] ]. 
Details for expressing the densities and currents in the 
cylindrical coordinate representation can be found in 
Ref . [43| • The pairing field is treated by using the density- 
dependent contact interaction (lij . 
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where Qo(r) denotes the isoscalar density of the ground 
state and P a the spin exchange operator. Assuming time- 
reversal symmetry and reflection symmetry with respect 
to the x~y plane, we have to solve for positive f2 and pos- 
itive z only, being the z— component of the angular mo- 
mentum j. We use the lattice mesh size Ap = Az = 0.6 
fm and a box boundary condition at p m ax = 9.9 fm, 
z max = 12 fm. The differential operators are repre- 
sented by use of the 11-point formula of Finite Difference 
Method. Because the parity and fl are good quantum 
numbers in the present calculation scheme, we have only 
to diagonalize the HFB Hamiltonian |T]) for each f2 7r sec- 
tor. The quasiparticle energy is cut off at E qPtCUt — 60 
MeV and the quasiparticle states up to O w = 15/2 ± are 
included. 

The pairing strength parameter t' is determined so 
as to reproduce the experimental pairing gap of 34 Mg 
(A eX p = 1-7 MeV) obtained by the three-point for- 
mula [H|. The strength t' = -295 MeV fm 3 for the 
mixed-type interaction (t' 3 — — 18.75^) [46| with 7=1 
leads to the pairing gap (A v ) = 1.71 MeV in 34 Mg. 

Using the quasiparticle basis obtained as a self- 
consistent solution of the HFB equations |T]), we solve 
the QRPA equation in the matrix formulation [53] 
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The residual interaction in the particle-hole (p-h) channel 
appearing in the QRPA matrices A and B is derived from 
the Skyrme density functional. We neglect the spin-orbit 
interaction term C^ J as well as the Coulomb interaction 
to reduce the computing time in the QRPA calculation. 
We also drop the so-called "J 2 " term Cj both in the HFB 
and QRPA calculations. The residual interaction in the 
particle-particle (p-p) channel is derived from the pairing 
functional constructed with the density-dependent con- 
tact interaction ([2]). 

Because the full self-consistency between the static 
mean-field calculation and the dynamical calculation is 
broken by the above neglected terms, we renormalize 
the residual interaction in the p-h channel by an over- 
all factor / p h to get the spurious mode. We cut the 
two-quasiparticle (2qp) space at E a + Ep < 60 MeV 
due to the excessively demanding computer memory size 
and computing time for the model space consistent with 
that adopted in the HFB calculation; 2_E qp cu t = 120 



TABLE I: Ground state properties of 36 > 38 ' 40 Mg obtained 
by the deformed HFB calculation with the SkM* interac- 
tion and the mixed-type pairing interaction. Chemical po- 
tentials, deformation parameters, average pairing gaps, root- 
mean-square radii for neutrons and protons are listed. The 
average pairing gaps of protons are zero in these isotopes. 
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-3.24 
-21.0 
0.31 
0.39 
1.71 
3.59 
3.18 



-2.41 

-23.7 

0.29 

0.38 

1.64 

3.67 

3.20 



-1.56 

-24.4 

0.28 

0.36 

1.49 

3.76 

3.22 



MeV. Accordingly, we need another factor / pp for the p- 
p channel. See Ref. (39| for details of determination of 
the normalization factors. In the present calculation, the 
dimension of the QRPA matrix ([3]) for the K n = 0~ ex- 
citation in 40 Mg is about 17 100, and the memory size is 
24.4 GB. The normalization factors are / p h = 1.06, and 
/pp = 1.21. 



III. RESULTS OF THE CALCULATION 

We summarize in Table Q] the ground state proper- 
ties. The neutron-rich Mg isotopes under investigation 
are prolately deformed. This is consistent with the re- 
sults calculated using the Skyrme SIII interaction [33j ] . 
The Gogny-HFB calculation using the D1S interaction 
suggested the shape coexistence in 38 > 40 Mg 36] . We 
can see that the neutron skin develops as approaching 
the drip line; the difference in neutron and proton radii 
y/ (r 2 )i, - x/ir 2 }^ = 0.41 fm in 36 Mg changes to 0.54 fm 
in 40 Mg. 

Figure [1] shows the response functions for the isovector 
(IV) dipole excitation in neutron-rich Mg isotopes. The 
IV dipole operator used in the present calculation is 



P ™k = e^r > :nY 1K (h) - e- A y'nY 1K {h), (4) 
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and the response function is calculated as 
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The Giant Dipole Resonance (GDR) appearing at 15 — 
25 MeV shows a deformation splitting for the K v = 0~ 
and 1~ excitations. In the lower energy region, we can see 
a bump structure above the neutron-emission threshold 
energy. 

Figure [2] shows the partial sum of the energy weighted 
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FIG. 1: (Color online) Response functions for the isovector (IV) dipole operator Q in 36,38 ' 40 Mg. The dotted, dashed and 
solid lines correspond to the K* = 0~, K* = 1~ and total responses, respectively. For the if = 1~ response, the transition 
strengths for the if = ±1~ states are summed up. The transition strengths are smeared by using F — 2 MeV. The arrows 
indicate the one-neutron continuum threshold E t h,in = |A| + minE a and the two-neutron continuum threshold £th,2n = 2|A|. 
In 40 Mg, these two continuum-threshold energies are almost degenerated. 
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FIG. 2: Energy weighted sum of the IV dipole strength func- 
tion in 40 Mg. The horizontal lines show the classical Thomas- 
Reiche-Kuhn (TRK), and the RPA sum rule values including 
the enhancement factor, mi = m\ (1 + k) (« = 0.33 in 40 Mg 
with the SkM* interaction ). 



strength defined as 

W{E X )= £&*l(*&|0>| 2 . (6) 
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For the IV dipole mode in °Mg, the calculated sum up to 
60 MeV reaches 98.5% of the EWSR value including the 
enhancement factor, mi = mf(l + k) where k = 0.33. 
The IV dipole strength below 10 MeV exhausts about 
6.0% of the classical Thomas-Reiche-Kuhn (TRK) sum 
rule. In 36 Mg and 38 Mg, the summed transition strength 
up to 10 MeV exhausts about 3.6% and 4.8% of the TRK 
sum rule, respectively. 

In what follows, the low-energy dipole excitations are 
investigated in detail. 



A. Low lying states in Mg 

Due to the deformation, the strength distribution and 
microscopic structure of the K w = 0~ and l - excitations 
are different. Figures [3Jc) and [3£f ) show the IV dipole 
strengths in the lower energy region for the K n = 0~ and 



TABLE II: QRPA amplitudes for the if = 1" state in 
40 Mg at 8.22 MeV. This mode has the isovector (IV) dipole 
strength B(IV1) = 6.24 x 10~Vfm 2 , the isoscalar (IS) oc- 
tupole strength B(IS3) = 1.44 x 10 3 fm 6 , and the sum of 
backward- going amplitude Y. Y lp = 5.56 x 10~ 3 . The single- 
quasiparticle levels are labeled with the asymptotic quantum 
numbers [JVn3A]fl. Only components with — Y^g > 0.02 
are listed. Two-quasiparticle excitation energies are given by 
E a + E/3 in MeV and two-quasiparticle dipole transition ma- 
trix elements D\ tCl p in e fm, the octupole transition matrix 
elements Oi ya /3 in fm 3 . In the rows (e) and (j), the label u3/2 + 
denotes a non-resonant discretized continuum state of neutron 
Q, w — 3/2 + level. The quasiparticle resonance of the hole-like 
^[200] 1/2 level is described by three discretized states in the 
present box size, and the level in the row (b) is the same state 
as in (g). 
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FIG. 3: IV dipole transition strengths in 36 ' 38 ' 40 Mg for the K* — CP (upper) and V (lower) excitations. 



1 excitations in 40 Mg. 

First, we are going to discuss the structure of the 
K n = 1~ excitations. At huj — 8.22 MeV, we can see a 
prominent peak possessing the large transition strength 
in Fig. H(f). We made a detailed analysis of this eigen- 
mode and show in Table |TT] its microscopic structure. 
This state is generated by many 2qp excitations. Among 
the 2qp excitations listed in Table HT1 the 2qp excitation 
of (d) and (h) is the particle-particle like excitation, and 
that of (g) is the hole-hole like excitation. These 2qp ex- 
citations never participate in generating the RPA mode 
in the absence of the pairing correlations. 

In order to understand the spatial structure of the state 
at 8.22 MeV, we show in the lower panel of Fig. 0] the 
transition density to this state. The transition density 
has the following features: Around the surface and the 
inside of the nucleus, the protons and neutrons oscillate 
with coherence. Outside of the nucleus, the neutrons only 
oscillate and the oscillation of the internal core and the 
neutron outside is in out-of-phase. This corresponds to 
the classical picture of an oscillation of the neutron skin 
against the isoscalar core along p— axis; the axis perpen- 
dicular to the symmetry axis. 

We are going to move on to the K v = 0~ excita- 
tions. Above the threshold energy, we can see several 
states possessing enhanced strengths in Fig. [3jc). The 
states at 4.83 MeV and at 6.20 MeV have large transi- 
tion strengths. The state at 4.83 MeV is generated dom- 
inantly by the 2qp excitation of i/[310] 1/2 ® i/ [440] 1/2 
(4.75 MeV for the 2qp excitation energy) with a weight, 
X 2 - Y 2 , of 0.88. The state at 6.20 MeV is gener- 
ated predominantly by the 2qp excitation of ^[310] 1/2 
and the discretized state of tl* = 1/2+ (6.26 MeV) 
with a weight of 0.75 and slightly by the 2qp excita- 
tions of i/[312]3/2 ® z/[411]3/2 (6.44 MeV) with 0.08 and 
^[200]l/2® z/[310]l/2 (7.22 MeV) with 0.03. 

We can see another prominent peak at 9.40 MeV. Mi- 
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FIG. 4: (Color online) Transition densities of protons and 
neutrons in 40 Mg to the K* = 0" state at 9.40 MeV (upper) 
and the K n = 1" state at 8.22 MeV (lower). Solid and dotted 
lines indicate positive and negative transition densities, and 
the contour lines are plotted at intervals of 3 x 10~ 4 fm -3 . The 
thick solid lines indicate the neutron and proton half density, 
0.055 fm~ 3 and 0.032 fin -3 , respectively. 



croscopic structure of this state is summarized in Ta- 
ble IIII1 This state is generated by many 2qp excita- 
tions as well as the K v = 1~ state at 8.22 MeV rep- 
resenting the pygmy dipole mode. The contribution of 
the qp excitation into the non-resonant continuum state 
is larger than the pygmy K v = 1~ state. The transi- 
tion density to this state is shown in the upper panel 
of Fig. [D Although the transition density represents the 
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TABLE III: Same as Table HQ but for the K w = 0~ state at 
9.40 MeV This mode has the IV dipole strength B(IV1) = 
7.99 x 10" 2 e 2 fm 2 , the IS octupole strength B(IS3) = 23.8 
fm 6 , and the sum of backward-going amplitude ^t»/3 = 
8.76 x 1CT 4 . In the rows (b), (d), (f) and (i), the label I/1/2+ 
and Z/3/2+ denote non-resonant discretized continuum states 
of neutron = l/2 + and 3/2 + levels. The quasiparticle 
resonance of the hole-like v [330] 1/2 level is described by two 
discretized states in the present box size, and the level in the 
row (g) is the same state as in (i). 
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FIG. 5: (Color online) Same as Fig. H but for the K n = 1" 
state at 8.67 MeV in 36 Mg (upper) and the state at 8.32 MeV 
in 38 Mg (lower). The neutron and proton half densities are 
0.051 fm -3 and 0.034 fm" 3 in 36 Mg and 0.053 fm* 3 and 0.033 



pygmy dipole character, an oscillation of the neutron skin 
against the isoscalar core along the symmetry axis, the 
amplitude is smaller than that of the transition density 
to the K* = 1" state at 8.22 MeV. 



B. Low lying states in 36 Mg and 38 Mg 

Figures EHa), E^b), OJd) and [3Je) show the transition 
strengths for the IV dipole excitation in the lower energy 
region in 36 Mg and 38 Mg. For the K n — 0~ excitations, 
we can see a peak at around 9 — 10 MeV both in 36 Mg and 
38 Mg. The state at 9.85 MeV in 36 Mg is generated by 
the superposition of many 2qp excitations; among them 
the 2qp excitations of i/[200]l/2 ® i/[330]l/2 (8.94 MeV) 
with a weight of 0.25, z/[330]l/2 ® z/[211]l/2 (10.3 MeV) 
with 0.17, and j/[202]5/2 <8> i/[312]5/2 (10.3 MeV) with 
0.13 have large contributions. The state at 9.69 MeV in 
38 Mg is also generated by many 2qp excitations. The 2qp 
excitations of i/[440]l/2(g>i/[321]l/2 (9.85 MeV) with 0.27, 
i/[330]l/2®i/[440]l/2 (8.83 MeV) with 0.15, and i/[200]® 
i/[330]l/2 (10.4 MeV) with 0.10 have large contributions. 
The transition densities to these states have a similar 
spatial structure to the transition density to the IC = 0~ 
state at 9.40 MeV in 40 Mg. 

For the — 1~ excitation, we can see a peak at 
around 8-9 MeV both in 36 Mg and 38 Mg. The state 
at 8.67 MeV in 36 Mg is generated mainly by the 2qp 
excitations of i/[202]3/2 <8 t/[321]l/2 (8.69 MeV) with a 
weight of 0.40 and i/[211]1/2®i/[310]1/2 (8.50 MeV) with 
0.37. The state at 8.32 MeV in 38 Mg is generated pre- 
dominantly by the 2qp excitation of z/[202]3/2(g)i/[321]l/2 



(8.33 MeV) with a weight of 0.69. 

Figure [5] shows the transition densities to the = 1~ 
states in 36 Mg and 38 Mg. These states have a structure 
of the neutron-skin oscillation against the isoscalar core 
similarly to the pygmy state in 40 Mg. As approaching 
the drip line, the neutron transition density has more spa- 
tially extended structure, and the amplitude is also larger 
whereas the proton transition density is unchanged. 



IV. COUPLING AMONG DIFFERENT MODES 
OF EXCITATION 

In a deformed system, the angular momentum is no 
longer a good quantum number. Therefore, we could ex- 
pect the coupling between the dipole and octupole modes 
of excitation in 40 Mg as one of the unique features of the 
negative-parity excitation modes in a deformed system. 
Figure [H^d) shows the K v = 1~ isoscalar (IS) octupole 
transition strengths. The lowest state at 4.68 MeV and 
the state at 5.79 MeV have enhanced octupole transi- 
tion strengths of 4980 fm 6 and 4990 fm 6 . The state at 
4.68 MeV is generated dominantly by the 2qp excitation 
of (i) i/[310]l/2 ® z/[440]l/2 (4.75 MeV) with a weight 
of 0.83, and the state at 5.79 MeV is generated by the 
2qp excitations of (ii) i/[310]l/2 <g> i/[411]3/2 (5.90 MeV) 
with 0.61 and (iii) i/[301]l/2<g>i/[440]l/2 (5.63 MeV) with 
0.11. The state at 4.68 MeV has a similar structure to 
the K v = 0~ state at 4.83 MeV. The excitation ener- 
gies do not change so much with respect to the unper- 
turbed 2qp excitation energies. Nevertheless, the transi- 
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FIG. 6: K 71 = 0~ (left) and 1~ (right) transition strengths for 
the IV dipole (upper), isoscalar (IS) octupole (middle) and IS 
dipole (lower) excitations in 40 Mg. 



tion strengths become large. This is because the unper- 
turbed transition strengths of the 2qp excitations of (i) 
and (ii) are quite large, (i) 858 fm 6 and (ii) 2070 fm 6 , as a 
consequence of the spatial extension of the quasiparticle 
wave functions around the Fermi level. 

We can see an appreciable coupling between the dipole 
and octupole excitations also for the K* = 1" pygmy 
state at 8.22 MeV possessing the IS octupole transition 
strength of 1440 fm 6 . 

Figure [7] shows the partial sum of the energy weighted 
strength for the IS octupole excitations. The spurious 
component of the center-of-mass motion is subtracted 
for the EWSR values for the K n = 0~ and 1~ excita- 
tions. (See Appendix for the effect of the cm. motion on 
the IS octupole excitations in a deformed system.) The 
summed octupole transition strengths up to 10 MeV ex- 
hausts about 26.6%(27.3%) of the EWSR value for the IS 
octupole A' 71 " = Q _ (l~) excitation. The individual eigen- 
states obtained in the present calculation scheme do not 
represent the collective nature. However, concentration 
of the transition strengths in the low energy region would 
be one of the unique features in drip-line nuclei. Inves- 
tigation of this unique feature is challenging in a more 
sophisticated framework that is able to handle the cou- 
pling to the continuum in a better way. 
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FIG. 7: (Color online) Same as Fig. [2]but for the IS octupole 
excitations. 



Enhancement of the transition strengths in the lower 
energy region for the IS compressional dipole excitation 
in unstable nuclei was also predicted in Ref. [48[. Since 
the pygmy dipole mode discussed above has both an IS 
and an IV structure, we can expect the enhancement of 
the strength for the IS dipole excitation as well as for the 
IV dipole excitation in Mg isotopes. 

Figure [8] shows the response functions for the IS com- 
pressional dipole mode. The IS dipole operator used in 
the present calculation is 



where 



pis _ 



V = 



E 




(A = 0) 
(A = ±l). 



(7) 



(8) 



The correction factor r\ originally discussed for a spherical 
system to subtract the spurious component of the cm. 
motion [49| was extended for a deformed system, and 
coincides with ry = | (r 2 ) in the spherical limit [39] . 

We can see a tremendous enhancement of the transi- 
tion strengths above the neutron-emission threshold en- 
ergy up to ~ 15 MeV, where the IV dipole strengths start 
to have a resonance peak. The enhancement of the tran- 
sition strengths are found not only for the pygmy modes 
but also for the less collective low-lying states possessing 
large IS octupole transition strengths as shown in Fig. [5] 

We can see a significant coupling among the IV dipole, 
IS octupole and IS compressional dipole modes in de- 
formed Mg isotopes close to the drip line. It is thus quite 
interesting to investigate systematically in a wide mass 
region of nuclei, where we can find the nuclear deforma- 
tion, not only appearance of the pygmy mode but also 
coupling between the dipole and octupole excitations and 
enhancement of the transition strengths for the IS com- 
pressional dipole and IS octupole excitations in the low 
energy region. In neutron-rich nuclei close to the drip 
line, the low lying modes are embedded into the contin- 
uum as shown in Fig. [SI Therefore, it is strongly desirable 
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E (MeV) 

FIG. 8: (Color online) Same as Fig. [1] but for the IS compressional dipole operator (0. 



to develop the microscopic framework that enables us 
to investigate the continuum effects on excitation modes 
and collectivity of the low lying resonance in a quantita- 
tive manner. Developing the QRPA on top of the HFB 
in Gamow basis [5011 is one of the ultimate ways and has 
been undertaken [5l| . The present work is considered to 
be a step toward the long-range plan of the microscopic 
and quantitative description of the collective modes of 
excitation in deformed drip-line nuclei, and the present 
results show that it is worthwhile to investigate the es- 
caping widths of these states and collectivity of the low 
lying resonances in a sophisticated framework. 
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V. SUMMARY 

We made a detailed analysis of the low-lying dipole 
states in deformed neutron-rich Mg isotopes close to the 
drip line, 36 Mg, 38 Mg and 40 Mg, by using the deformed 
QRPA employing the Skyrme SkM* and the local mixed- 
type pairing energy-density functionals. 

Above the neutron-emission threshold, we obtained 
several eigenstates having large transition strengths for 
the IV dipole excitation. Because of the deformation, 
excitation modes with different K quantum numbers 
have different excitation energies and microscopic struc- 
tures. We obtained the collective pygmy-dipole modes at 
around 9 — 10 MeV for the K v = - excitation and at 
around 8 — 9 MeV for the K w = 1~ excitation. These 
pygmy dipole modes are generated by many 2qp excita- 
tions. 

Since the pygmy dipole modes obtained in Mg isotopes 
have both an IS and an IV structure, it has been found 
that the transition strengths for the IS compressional 
dipole excitation are tremendously enhanced in the lower 
energy region. Furthermore, because of the mixing of dif- 
ferent angular momenta in a deformed system, we found 
significant coupling among the IV dipole, IS octupole and 
IS compressional dipole modes in the low energy region. 



APPENDIX: SUBTRACTION OF THE 
CENTER-OF-MASS MOTION FROM THE 
ISOSCALAR OCTUPOLE OPERATOR 

Because of the deformation, the = 0~ and 1~ 
octupole excitations contain the spurious center-of-mass 
motion. We deal with this problem by using an operator 

M 3Q {r) = r 3 Y 3Q - V rY w (A.l) 

= \[^~{^ - 3zp 2 - 7]'z), (A.2) 

where rj = y/Yljlr\ for the K n = 0~ channel by follow- 
ing the discussion in Ref. [49( • The correction factor rf is 
determined by the condition of the translational invari- 
ance. 

The vibrating density associated with the external field 
V(r,t) = a(t)M 30 (r) (A.3) 
can be expressed to first order in a(t) by 
«e(r > t)=aV-(eoVM 3 o) (A.4) 

= f£ a { { - 6zp) Y P + (6z2 - 7,1 - 3p2) lL} g °- 

(A.5) 

Here Qq is the ground-state density, and the variables t 
and r are omitted for simplicity. 
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The condition of the translational invariance of the sys- 
tem 



gives 



J dr6erYf = 



(A.6) 



(A.7) 



The similar procedure is taken for the K~* = 1 chan- 
nel by using an operator 



M 3 i = r 3 Y 31 - W Y tl 



21 



(4z 2 p-p 3 -r/p)e i +, 



64tt 

where rf = y/8/7rj. It gives 

rf K=1 =4(z 2 )-2{p 2 ). 



(A.8) 
(A.9) 

(A.10) 



In the spherical limit, the correction factors both for 
the K n = CP and 1~ excitations vanish. This is rea- 
sonable because the octupole excitations decouple to the 
spurious cm. motion in a spherical system. 

The EWSR values for the IS octupole excitations 



(A.11) 



are given by 

EWSR(A = 3, K = 0) 



(A.12) 



EWSR(A = 3, K = 1) 

= ^ AX ^t (16( " 4) + Hp4) + 16<zV » ' (A ' 13) 
EWSR(A = 3, K = 2) 

h 2 . 105 



A x ^ (8(p 2 z 2 ) + (p 4 )) 



2m 32tt 
EWSR(A = 3, K = 3) 

ft 2 315 4 
~ 2^ X 32^ P ' 



(A.14) 



(A.15) 



The EWSR values for the K 71 = and 1 excitations 
are corrected as 



EWSR cor (A = 3,K = 0) 

= EWSR(A = 3, K = o; 
EWSR cor (A = 3,if = 1) 



= EWSR(A = 3, K = 0) - £-A x j^Vk=o, (A.16) 



EWSR(A = 3, Jf = 1) - |-A x ^-^=i (A.17) 



by subtracting the spurious component of the cm. mo- 
tion. 
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